A gauge model for right handed neutrinos as dark matter 
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B ~ L is known to be a symmetry somewhat linked to the origin of neutrino masses, which turns 
out to be anomaly free upon the sole introduction of right handed neutrinos. We suggest a simple 
extension of the electroweak group, SU{2)l x U{1)y x U{1)b-l, where the breaking of U{1)b-l 
symmetry provides masses for right handed neutrinos at an acceptable range for them to be Dark 
Matter candidates. We review the cosmological constraints for this type of models, in order to find 
the restrictions on the parameters of the model. We study the contributions due to the B — L 
current in the process e"'"e~ — >■ e"'"e~ at low energies. 

PACS numbers: 14.60.St, 95.35.+d 



I. INTRODUCTION 

There is compelling evidence that the Standard Model 
(SM) of strong and electrovifeak interactions is not com- 
plete. There are several experimental and observational 
facts that cannot be explained within the SM. These 
are neutrino oscillations [iL the presence of Dark Mat- 
ter (DM) in the Universe 3|, the baryon asymmetry of 
the Universe, its flatness, and the existence of cosmolog- 
ical perturbations necessary for structure formation Q. 
Neutrino oscillations give indirect proof that neutrinos 
are massive. Their mass, however, needs that the SM 
be extended at least to incorporate right handed neutri- 
nos (RHN), and the simplest way to understand the tiny 
scales required for neutrino oscillations is to implement 
the see-saw mechanism 0, [H| . The simplest extension to 
the SM that includes the seesaw Lagrangian structure is 
written as 



CuMSM = CsM + NIl^^,-1^'NI 



N'lNi + h.c, (1) 



where Csm is the Lagrangian of the Standard Model, 
$i — ey $* and L^, (a = e,iJ,,T) are the Higgs and lep- 
ton doublets, respectively, and both Dirac {M^ = ^{^)) 
and Majorana (M/) masses for neutrinos are introduced, 
where / = 1,2,3 stands for the number of RHN species. 
The theory based on this Lagrangian has been called the 
j/-Minimal Standard Model (i^MSM) In comparison 
with the SM, the j/MSM contains 18 new parameters and 
interestingly enough, it appears to have the potential to 
provide some explanations to the above mentioned prob- 
lems 0. 

The new parameters of the i^MSM can describe any 
pattern of masses and mixings of active neutrinos, which 
is characterized by 9 parameters only. Despite of this 
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freedom, the absolute scale of active neutrino masses can 
be pin pointed in the j/MSM from cosmology and astro- 
physics of dark matter particles [6l-[l0j : one of the active 
neutrinos must have a mass smaller than 0(10"^) eV. 
The choice of the small mass scale for singlet fermions 
leads to small values of the Yukawa coupling constants, on 
the level of 10^^ — 10^^^, which is crucial for the ex- 
planation of dark matter and baryon asymmetry of the 
Universe. 

Although the i/MSM does not have any extra stable 
particle in comparison with the SM, the lightest singlet 
fermion, A''! , may have a life-time greatly exceeding the 
age of the Universe and thus play a role of a dark matter 
particle ^MB- 

DM sterile neutrinos can be produced 
in the early Universe via active-sterile neutrino transi- 
tion (llj : via resonant active-sterile neutrino oscillations 
in the presence of lepton asymmetries ; or during in- 
flation. DM sterile neutrino may also have other inter- 
esting cosmological applications [16| , in particular for the 
understanding of baryon asymmetry of the Universe. On 
the other hand, the baryon (B) and lepton (L) numbers 
are not conserved in the z^MSM. The lepton number is 
violated by the Majorana neutrino masses, while B + L 
is broken by its anomaly. As a result, the sphaleron pro- 
cesses with baryon number non-conservation [l7j are in 
thermal equilibrium for temperatures 100 GeV < T < 
1012 GeV. As for CP-breaking, the i/MSM contains 6 CP- 
violating phases in the lepton sector and a Kobayashi- 
Maskawa phase in the quark sector. This makes two 
of the Sakharov conditions [l^ for baryogenesis satis- 
fied. Similarly to the SM, this theory does not have 
an electroweak phase transition with allowed values for 
the Higgs mass [l^, making impossible the electroweak 
baryogenesis, associated with the non-equilibrium bubble 
expansion. However, the i/MSM contains extra degrees 
of freedom -sterile neutrinos- which may be out of ther- 
mal equilibrium exactly because their Yukawa couplings 
to ordinary fermions are very small. The latter fact is 
a key point for the baryogenesis in the :/MSM [2^ [2l| , 
ensu ring the validity of the third Sakharov condition. In 
Ref. [21| it was shown that the i/MSM can provide simul- 
taneous solution to the problem of neutrino oscillations. 
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dark matter and baryon asymmetry of the Universe. 

However, this theory breaks exphcitly the global sym- 
metry B — L, due to the term NfNj, which also rep- 
resents the introduction by hand of a fermion mass in 
the scheme, against the established lore in the standard 
theory, which claims that fermion as well as gauge boson 
masses are due to the Higgs mechanism. We consider this 
as a serious theoretical draw back of the model. This is, of 
course, a known fact that usually drives theorist to con- 
sider the left-right extensions to the standard model (see 
for instance Ref. [22*1 ). Nevertheless, which seems quite 
appealing to exploit is the fact that solely introducing 
three families of RHN's in the SM particle content comes 
with the extra free ingredient of making B—L an anomaly 
free symmetry, without any further requirement. Inter- 
estingly enough, keeping the physical connection between 
the mass of RHN's and the Higgs Mechanism, suggests 
the simplest gauge group SU{2)l x U{1)y x U{1)b-l as 
a natural symmetry which extends the electroweak group 
of SM, and which should contains most of the features 
of the j/MSM. In here, B — L turns out to be unrelated 
to hypercharge, and thus, the model offers a completely 
different route to understand the origin of see-saw terms, 
and, of course, the gauge origin of the ^MSM, than that 
provided by left-right type models. To explore some of 
the features of such gauge extension is the main goal of 
the present work. 

This paper is organized as follows, first we present the 
model SU{2)l x U{1)y x U{1)b-l and give the justifi- 
cations for adding the U{1)b-l gauge symmetry to the 
electroweak group. Then, as the new interactions intro- 
duced on the model for RHN's do require us to analyze 
their potential over the former i^MSM considerations for 
DM, we check the cosmological constraints on the model. 
And finally, we discuss the considerations under which 
the model could be consistent to collider results. In par- 
ticular we focus on low energy Bhabha scattering, which 
could constrain a weakly interacting light B — L boson, 
imposing preliminary bounds on the gauge coupling to 
mass ratio. 



II. THE B-L MODEL 



The model under consideration is based in the gauge 
symmetry group SU{2)l x U{1)y x U{1)b-l, and, as we 
have already stated, it is motivated on the idea that the 
Higgs mechanism should be regarded as the fundamental 
mechanism for the generation of masses. Matter content 
is the same as for the SM, but for the addition of a RHN, 
TV, per family, and a SM singlet scalar, a, that is going to 
be used to implement spontaneous symmetry breaking. 



So, this model contains the following particles, 

L ^(2,-1,-1); 
CR ^ (1,-2,-1); 

N ^(1,0,-1); (2) 
H ^(2,1,0); 
a ^(1,0,2); 

where the numbers label the Weyl representation of 
{SU{2)l, U{1)y,U{1)b^l), and family indices arc under- 
stood. It is not difficult to check that, with this fermion 
content, B ~ L is anomaly free [23], and that Yukawa 
couplings, 

yiaN^Ni, (3) 

will generate the RHN mass terms of the i^MSM, M//2 — 
yi{<j), upon spontaneous symmetry breaking, thus, re- 
lating the RHN mass scale to the physical scale at which 
B—Lis being broken. It is important to notice that B—L 
in this model appears as a symmetry that is orthogonal 
to the hypercharge, and therefore with no relation, what- 
soever, to the electric charge. This is a distinctive feature 
that makes the study of the model worthy. 

Several models exist on the literature that had added 
to the SM, a U{1) gauge symmetry 24-30]; those models 
assume that the range of the breaking scale of the sym- 
metry is bigger than few hundred GeV ^3i|, H^l , which is 
usually due to the mixings among the associated B — L 
gauge bosons and the standard Z boson, for which the 
usual bounds on exotic Z' searches apply |33|]. Without 
the tree level mixing, it is clear those bounds do not apply 
straightforwardly. In contrast, the search oi B — L gauge 
boson, Zb~l, through dilepton, or dijet production in 
colliders, seems more promising. A complete analysis 
for extra Z' gauge bosons from a collider point of view 
can be found in Ref. ,34j. In particular, for a vector-like 
U{1)b~xL gauge boson coupled to left and right handed 
electrons, one can extract the lower bound due to LEP 
experiments, 

Mz' > \x\gz' X (6 TeV). (4) 

Thus, in a _B — i model there is still a wide range al- 
lowed for the mass of the extra gauge boson depending 
on the coupling. So, we can stretch the LEP limit on 
the ignorance of the value of gauge coupling associated 
to U{1)b-l- Therefore, the possibility that this gauge 
group might be broken at scales similar to electroweak 
scale is yet present, and worth of consideration. 

In the model described here, it is supposed the ex- 
istence of two Higgs fields, one which breaks the elec- 
troweak sector {SU{2)l x U{1)y), H{x), and the second 
which breaks the U{1)b-l symmetry, a{x). The most 
general lagrangian which should be used to implement 
symmetry breaking is 

C = {d,H)Hd'^H)~mj,\H\' + ^\H\'' + 

id,a)Hd^a) - ml\a\^ + + (5) 



where a, /3 and S are constant parameters. The parame- 
ter S is related with the mixing between the electroweak 
and the B — L groups and we are going to work in the 
hmit where 5 <^ p. This is the most simple realization 
of U{1)b-Li where the electroweak gauge group is not 
mixed with the B ~ L one, thus, the results presented in 
here would be considered as a first approximation of the 
full SU{2)l X U{1)y X U{1)b~l gauge theory. 

We still have not identified the breaking scale of 
U(1)b-l and, in this context, we will not be able to 
determine it, unless an experimental evidence of a S — L 
gauge boson appears. Hence, we have two possibilities for 
the breaking scale: i) if {H) <^ {a) we end up in the class 
of models where mzB-L ^ and the physics related to 
these gauge groups is decoupled; and ii) if {H) ^ (a), the 
physics involved in this scenario is such that the gauge 
bosons are not decoupled. Moreover, the non evidence 
of new physics hits also to this class of models, meaning 
that the probable effects might be suppressed in some 
manner. In this work, we consider that mzs-L ^ '^^Zj 
the motivation of this proposal is based on the analysis 
made by Ref. (ssl . IsgI ] , where it was pointed out that a 
RHN mass of the order of keV is a good candidate to be 
the DM particle, and, in this context, our above assump- 
tion might lead us to this DM scheme. 

By imposing the spontaneous breaking of the local 
U(1)b-l symmetry, a massive neutral gauge boson, 
Zb-Li emerges with a mass 



(6) 



where gB-L is the coupling constant associated to 
U{1)b-l and (cr) is the vacuum expectation value of the 
a field. In addition, Zb-l couples to matter through the 
coupling —gB-LZ^~^ Jg_i^, where the vector-like B — L 
current is given as. 



Jb-l = ^lY'^l + CLl^eL + eRY'eR, + N-yi^N . (7) 

With this simple structure, we can implement several 
searches oi B — L currents in different levels. First, a 
realistic model for RHN as dark matter must satisfy the 
cosmological limits which are set due to several obser- 
vations. So next, we are going to review the applicable 
cosmological constraints for this model. 



III. COSMOLOGICAL DARK MATTER 
CONSTRAINTS 

Let us start the discussion of the model by looking 
at the cosmological constraints. As it is pointed out in 
Ref. [stI, models with RHN as DM must satisfy the con- 
strains given mainly by: a) the structure formation re- 
quirement, b) X/7-ray observations and, c) the abun- 
dance of these particles at present time. 

If the mass of the DM candidate is too light, the ob- 
served structure of the Universe would be eliminated by 
a too hot DM. The most restrictive bound on the mass 




FIG. 1: 

decay. 




Feynman diagrams which contribute to A'^ 
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of the DM candidate is coming from the Lyman-a forest, 
where the bound constraints the velocity distribution of 
the DM particles from the effect of their free streaming on 
the formation of the structure on Lyman-a scales. The 
bound on our lightest RHN mass is thus given by [37| . 



Ml > 1.6 kcV 



(8) 



Therefore, we are going to work with this value, not to 
spoil the large structure formation. 

On the other hand, a general feature of this model is 
that the RHN, the Dark Matter candidate, does decay 
into SM particles. This characteristic appears basically 
due to the fact that RHN could transform into an active 
one via the mixing parameter. 



E 



Ml 



(9) 



where v is the SM Higgs vacuum expectation value. 
Thus, one can compute the decay rate of RHN into 
charged leptons plus a SM neutrino {£~^i~i') or into three 
SM neutrinos in order to check a natural DM require- 
ment: the lightest RHN should be long lived and, fur- 
thermore, its lifetime has to be longer than the age of 
the universe. Nevertheless, this constraint is less dom- 
inant compared to what it results from the X/7— rays 
observation. The radiative decay N — > 1/7, induced at 
one loop level (FIG. 1), produces a narrow line in the X- 
ray spectrum of astrophysical objects [ll,|3l|. The width 
of this decay is given by [39| . 



■ N- 



9aGl . 2 



10247r4 



sin^(26ii)Mf, 



(10) 



where a is the fine-structure constant, and Gp = 1.166 x 
10~^ GeV~^ is the Fermi constant. This width can re- 
strict the mixing parameter when we consider the upper 
bound found in Ref. [st'], which reads as. 



F < 9.9 X 10^27 ggp-i 



and it is translated to 



9l < 1.7 X 10" 



1.6 keVy 
Ml ) 



(11) 



(12) 



Within the context of the parameters of the model, this 
result can be seen as a bound to the Yukawa parameters. 
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Therefore, no restriction on the B — L couphng and the 
mass of the associated gauge boson can be made so far. 

Finally, the abundance of RHN at the present time can 
be calculated by using [37l |. 



n 



N 



135C(3) Ml 



So 



DM 



(13) 



where VLdm = 0.105 h-'^, so = 2889.2 cm'^, = 
1.05368 X 10~^ h'^ GeV cm^^, 5 the entropy release and 
g*/ is the effective number of degrees of freedom imme- 
diately after freeze-out . Therefore, for a 1.6 keV RHN 
mass, we have. 



n 



N 



DM 



10.75 

9*f 



Ml 



1.6 keV 



X 160, 



(14) 



where we have chosen g*/ = 10.75, this value corresponds 
to the case where there are only SM particles and the 
freeze-out have happened below 100 MeV. In order to 
not overdose the universe we need a dilution factor of 
the order of. 



S"- 160 X 



10.75 

fJ*f 



Ml 



1.6 keV 



(15) 



we can achieve this requirement by considering that heav- 
ier RHNs could generate such amount of entropy. Let us 
consider that only the heaviest right handed neutrino, 
iVa , is the author of such entropy. In the case where the 
entropy generation is large, the calculation gives [37j , 



0.76 X 



9n 



-1/4 

9* ruH 



2 g^^yVMpi 



(16) 



where = 2 corresponds to the number of degrees of 
freedom of N^, and g^ are the properly averaged ef- 
fective number of degrees of freedom during decay 
and at freeze-out, and mn is the heavy neutrino mass. 
Therefore, Eq. (IT5|) and p6)) can be solved in order 
to constraint m//, nevertheless, the entropy generation 
should end before the big bang nucleosynthesis occurs 
(BBN), saying, the reheating temperature should satisfy 
the limits [1351 0.7 < < 4 MeV. This temperature 
is approximatively [isj . 



27r^g* 
45 



^1/4 



VTM 



PI- 



(17) 



Combining Eq. (fTS]). ([T6| and (|17p . we find a constraint 
on rnu to be. 



rriH > 



Ml 



1.6 keV 



(2.5^ 15) GeV 



(18) 



Although, it has been pointed out by some authors 'st'] , 
that once we have a constraint on the mass of the heavy 
RHN, it is possible to get a bound on the mass of the 
extra gauge boson. In order to do so, one needs to make 



use of the following fact: whether a heavy RHN is decou- 
pled, the decoupling must happens when the freeze-out 
temperature is higher than the mass of this heavy RHN. 
So, we just need to estimate the freeze-out temperature 
by considering that the mean free path of RHN should 
be equal to the Hubble scale. In order to estimate the 
moment when RHN decouples to the SM particles, we 
can use that the processes involved {NiNi <->■ e~^e~ , etc.) 
which are analogous to those for the usual SM neutrinos; 
the only difference appears in the coupling and in the 
mass of the B — L gauge boson, therefore, 



'NiNi 



W 



gs-L M, 
g M 



(19) 



where g is the SU{2)l gauge coupling, and AIw is the 
mass of the W boson. And so, the freeze-out tempera- 
ture, T/, would be 

T,^gy^(^^Y\l^2)MeY. (20) 

\gB-L Mw J 

Hence, one gets a bound for B — L gauge boson mass 
which reads as, 



M > 



9 



(6 10) TeV, 



(21) 



which is remarkably close to the limit set by collider 
searches by a few percentage. Furthermore, this re- 
sult matches in the limit when the gauge bosons cou- 
ple equally, to the result presented in Ref. [s^l- Nev- 
ertheless, we have no idea about the absolute scale of 
9b-l so far. This is important since it defines whether 
B — L physics could be relevant at any given scale. For 
instance, for a sufficiently small value of the coupling, say 
gs-L ^ 0(10"^), one could resize the lower bounds on 
M down to few GeV scale, suggesting that former exper- 
iments should have be sensible to the involved physics. 
It is, therefore, of a clear interesting to know how big the 
value oi gs-L is, in order to know whether its effects can 
be seen in future experiments. 



IV. CONSTRAINING THE MODEL USING 
PETRA EXPERIMENT 



Past low energy experiments can help us to upper 
bounding the value oi gp-L- A very clean channel where 
an extra neutral gauge boson can appear is in e+e~ scat- 
tering, and this reaction was studied by the PETRA ex- 
periment. We now reanalyze PETRA results in order to 
explore for a possible low scale mass of Zb-l- This limit 
would give us a hint about the values of gs-L or (a), and 
therefore we can stretch the limit in order to fulfill the 
DM constraints. Although, LEP and Tevatron experi- 
ments, and eventually LHC data, can also be included in 
the analysis, the known bound already stated in Eq. 



5 




FIG. 2: Feynman diagrams for e — >■ e in the SU (2) _l x 
)7(l)y X U{1)b-l theory. 

get particularly soften for small values of the B — L cou- 
pling, as we mentioned above. Thus, we would expect 
little improvement of the combined analysis, whereas it 
would become more involved. For simplicity, we left this 



analysis for a further work. In the following, we will use 
this experiment to restrict whether the mass of the gauge 
boson or the gauge coupling. 

PETRA results had been well explained theoretically 
by using the simple Bhabha scattering process. The exis- 
tence of a light and weakly coupled Z^~^ can be strongly 
restricted by checking the contribution of whatever extra 
neutral gauge boson to the Bhabha cross section. So, we 
perform the computation of e~e^ — > e~ scattering to 
establish a bound on gs-L and M . The relevant Feyn- 
man diagrams for the process are pictured in FIG. 2. 

Thus, in the limit E ^ m, Bhabha's expression re- 
ceives a contribution from Zb-l and the result can be 
read as, 



da 

'dn 



da 

dn 



O 



9%-L 



(22) 



where (^)^ is the Bhabha's differential cross section and the correction due to the extra gauge boson is enclosed in 
the S function, such that. 



2(a;,y) 



-16(1 - v){-2\Q -f 4x(9 + 3y - - iy^) + Sy'^{-IA + 8y -f y'^)) 
[x - 4){x + 2 - 22/)(-216 + 8y2(-l4 + %y + y^)) ' 



(23) 




FIG. 3: Correction to PETRA experiments. On the left panel we plotted the results of the SU{2)l x U{1)y x U{1)b-l with 
the experimental data. On the right panel we plotted the ratio of the experimental measure and the theory described here, for 
-^s = 34.8 GeV, due to the experimental data is more accurate at this center of mass energy. 



A numerical analysis of (|22]) shows that for any value 
gB-L < 10"^ m El, PETRA resuhs can be satisfied 
independently of the value for M, provided M ^ E. 
Indeed, by using gB-L = 10^^, for instance, it is possible 
to reproduce the PETRA results [1^, as it is shown in 
the FIG. 3, for M = 90 GeV. 



less restrictive, meaning that we can allow masses of the 
Z' as low as 90 GeV and remain still consistent with 
the previous observations. If we consider, naively, that 
M — 90 GeV, we are also considering that the breaking 
scale oi U{1)b-l is occurring at 



For such an small value for gB-L, LEP limit becomes 



{a} - 3 TeV, 



(24) 



6 



which is a reachable scale for the LHC. 

On the other hand, the mass of the right handed neu- 
trino is not fully controlled by (cr) . Considering the Ma- 
jorana term in Eq. ([3]), we can still have a RHN mass of 
the order of keV by choosing the Yukawa parameter to 
be 2/1 2.5 X 10"^". For the heaviest RHN, a Yukawa 
larger than 10^'^ would be enough to satisfy the bound on 
Eq. (jl8p . Certainly, on this path, one gives up natural- 
ness on the see-saw mechanism, but gain the possibility 
of accounting for DM. 

V. CONCLUDING REMARKS 

In this paper we have analyzed a simple extension 
of the electroweak gauge group by adding an extra 
U{1)b-l- This combination of barion and lepton num- 
bers makes the model anomaly free. Moreover, with this 
combination of quantum numbers, we have been allowed 
to write the RHN field and, by choosing Yb-l = 2 for 
an extra Higgs field, cr, we have been able to write a 
Majorana mass term for neutrinos when this extra Higgs 
field acquires a vacuum expectation value. Spontaneous 
breaking of U{1)b-l gauge leaves a a massive Zb-l 
gauge boson, which related to the breaking scale, and the 
coupling constant associated to U{1)b-l- By adding the 
corresponding RHN fields, we can address the problem of 
the DM content in the universe according to Ref. (35l.[36|. 
Thus, we studied the minimum conditions to recover the 
conclusions done for z^MSM in the DM point of view. In 
order to say something about neutrinos as dark matter, 
we have analyzed the cosmological constraints. These 
constraints had required to the lightest RHN to have a 
mass of about few keV in order to satisfy the observed 



structure formation. The X/7— ray observation has set 
a bound on the mixing between active and sterile neu- 
trinos to be 61 < lO^*'. In order to produce the correct 
abundance of RHN as DM, we have restricted the mass 
of the heavy RHN to be of the order of few GeV. The 
imperative decoupling of the heavy RHN from the light- 
est one has set a limit on the mass of the Zb-l, which 
is comparable to that obtained from LEP data; never- 
theless, the relation found is not conclusive since we still 
have the possibility of a small B — L gauge coupling, such 
that, we could still have a consistent light Zb~l gauge 
boson. Thus, in order to find an upper bound to gB-L we 
analyzed the possibility that the net effect of U{1)b-l 
could be suppressed in the process e+e^ — > e+e^; in par- 
ticular, since it must appears as a contribution into the 
Bhabha scattering. The simplest experiment which we 
can compare to is PETRA which has already measured 
exhaustively this scattering process. A combined analy- 
sis for LEP, Tevatron and eventually LHC experiments 
is left for a future work. The result indicates that, if the 
couphng constant associated to U{1)b-l is of the order 
of 10~^ or less, the contribution given by Zb-l is indeed 
subleading to the level of being mostly negligible. With 
this condition, we are also imposing that the breaking 
of U{1)b-l should occurs at about 3 TeV or higher, a 
sizable value that come into the energy region of interest 
for the LHC. 
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